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ABSTRACT 

In this paper we make an investigation of the receptivity of boundary layer flows to Gbrtler 
vortex modes. A study by Denier, Hall & Seddougni (1991) of the generation of vortices 
by wall roughness elements concluded that such elements are extremely poor as mechanisms 
to stimulate short wavelength modes. That work also examined the equivalent problem 
pertaining to 0(1) wavelength modes but that analysis was in error. We re-examine this 
problem here and demonstrate how the form of the wall roughness is crucial in determining 
t le vortex stability characteristics downstream of the roughness. In particular we investigate 
■e cases of both isolated and distributed forcing functions and show that in general a 
distributed function is much more important in generating vortices than are either isolated 
roughness or free-stream disturbances. 
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1. Introduction 

Until quite recently little was known about the important problem of the receptivity of 
boundary layers to Gbrtler vortices. Perhaps the first work in this direction was performed 
by Hall (1990) who was concerned with the receptivity question in the case of vortices 
excited by free-stream disturbances. That paper also reviews the various stages of growth 
of Gortler vortices within boundary layers and the reader is referred to that article for an 
extended discussion of the linear and nonlinear characteristics of Gortler modes. 

The receptivity work of Hall (1990) used a linearised analysis and this was extended by 
Denier, Hall & Seddougui (1991) (hereafter referred to as DHS) in order to account for the 
processes by which Gortler vortices may be stimulated by wall roughness elements. DHS 
considered both isolated and distributed roughness patches and they calculated coupling 
coefficients which indicate how the amplitude of the forcing function is related to the 
induced vortex size. For small vortex wavelength the coupling coefficient is exponentially 

small and so it was concluded that this type of mode is unlikely to be directly excited by 
wall roughness. 

This finding concerning the receptivity of boundary layers to high wavenumber vortices 
was deduced using an analytical technique in which the flow structure was described by 
WKB type expansions. However, in general a numerical method is required in order to 
investigate the properties of 0(1) wavenumber Gortler vortices generated by roughness. 
Laplace transform solutions of the governing equations were used by DHS in order to 
calculate the structure of the vortices at small distances downstream of the roughness 
region. This structure could then be used as initial conditions for the numerical marching 
of the relevant equations in order to obtain the flow properties at 0(1) distances from 
the roughness. These initial conditions assume the form of similarity solutions and the 
disturbance structure as predicted by DHS implies that close to the roughness the vortex 
is effectively confined to a very thin layer adjacent to the bounding surface of the flow. 

The type of analysis outlined in DHS was also used by Hall & Morris (1992) in the 
context of describing how vortices develop in the flow over a heated flat plate when these 
modes are triggered either by localised or distributed wall roughness or by non-uniform 
wall heating. Both these papers concluded that for 0(1) wavelength modes a vortex field 
grows wit hm a wedge at a finite distance downstream of the non-uniformity in the flow. 

Subsequent investigations of the receptivity of Gortler vortices have suggested that the 
work of DHS for 0(1) wavenumber modes is erroneous for two reasons. First, the inversion 
of the Laplace transform solutions derived in DHS are technically more involved than as 
was presented. This in turn has implications for the similarity solution structures and we 
show below that the required structure is not confined in the thin wall layer predicted by 
DHS. Therefore the numerical solutions obtained by DHS by integrating the full equations 
were initiated with incorrect disturbance profiles and so the aim of the current work is to 
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reconsider tl.is 0(1) wavenumber study. This allows us to illustrate how the analysis and 
conclusions of DHS need to be amended. First we concentrate on the problem with an 
isolated roughness element and then demonstrate how useful deductions concerning the 
distributed situation may be obtained. 

To begin we shall give a brief derivation of the pertinent equations, fuller details of 
which may be found in DHS. Consider the flow of a viscous fluid over a wall of variable 
curvature. If Uoo is a typical flow velocity far from the wall, v is the kinematic viscosity of 
the fluid and L is a typical lengthscale over which the wall curvature varies then we define 

a Reynolds number Re by 

Re = UooL/u , ( L1 ) 

and we suppose that with respect to Cartesian axes x*,y*,z* the wall is defined by 

y* = L f Re~?g(x*/L) + ARe~?f (x* /L,Re?z * /!<)] , (!- 2 ) 

where A is a small constant. Variables (x,y,z) are defined according to 

(x,y, z) = (x* , Re* y* , Re^ z )/L , (1.3a) 

and the corresponding velocity vector is supposed to be given by 

(u + ,u + ,ui + ) = (u* , Re% v* , Re% w*) /Uoo' O- 36 ) 


Attention is confined to the limit Re — ► oo and then 

p + = p(x) + Ai2e _ 2p(ar,y,z) + 0(A 2 ), (u + ,u + ,«> + ) = (u,v,0) + A(u,v,w) + O(A 2 ) 

(1.4a, b) 

where p + is the fluid pressure scaled on pU^ with p the fluid density. The basic flow 
velocities fi,» are functions of * and y alone whilst u,v,w,p depend on all three spatial 
co-ordinates. Expressions (1.4a, 6) are substituted into the linearised Navier-Stokes and 
continuity equations and the resulting forms need to be solved subject to no-slip conditions 
on the wall y = g and need also match to the free stream fluid velocity as y - oo. As 
in DHS we shall restrict our analysis to Blasius flow as a specific example although it is 
straightforward to generalise the work to other basic boundary layers. 

It is convenient to apply a suitable Prandtl transformation to the linearised disturbance 

equations to give 

i I + 5„+m, = 0, (1'5«) 


UU X + VUy + VUy + uu x — Uyy + 

UV X + VVy + UV X + VVy + Gx(x)uu = -Py + Vyy + V Z z, 


(1.56) 

(1.5c) 
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uw z + VU)y = -p z + w + 


w 


ZZ > 


“ = —fu y , V = W = 0 

i l,v,w —* 0 

whilst the basic flow quantities satisfy 


on 


y = o, 


as 


y -*• oo, 


(1.5rf) 

(1.5e) 

(1.5/) 


u x + Vy= 0, uu x + VUy = -p x + u yy , 


u — v = 0 on y = 0; 


u 


1 


as y — * oo. 


In equation (1.5c) the quantity 2(dV«fc 2 ) has been replaced by G X (x) where X and 
G are referred to as the wall curvature and Gdrtler number respectively. It is convenient 
to assume that the function f(x,z) (which determines the geometry of the roughness 
on the wall) may be written as / = f( x )q(z) and then equations (1.5) may be Fourier 
transformed in z. If denote the Fourier transforms of u,v,... etc. then it is 

convenient to eliminate w and p from the transformed versions of (1.5). If k denotes the 
transform variable then we obtain 


UU X + VUy + VUy + UU X — Uyy 


~ u yy ~ k 2 u, 


(1.6a) 


{ u xyy + k 4 + k 2 Vy} V + V x U yy + {u xyy + k 2 V x + k 2 X Gu} U 


+ {uyy - ud 2 /dy 2 + fc 2 u} V X + 2 { u xy + u x d/dy } u x (1.66) 

+ v yyyy ~ vv yyy - {u s + 2k 2 } v yy + {u xy -f- k 2 v } v y = 0, 

linking u and v whilst w and p may be then deduced from the transformed versions of 
(1.5a, d) 

u x +v y + ikw = 0, (1.6c) 


uw x + vw y = —ikp + U)yy — k' 


w. 


(1.6(f) 


should these quantities be desired. On defining F(x) = fi,(x,0)/(*) it is found that 
equations (1.6) need to be solved subject to 


u — q(k)F(x ), v = w = 0, 


on 


y = 0, 


u, v , w 


0 


as 


y 


oo. 


(1.6e) 

( 1 - 6 /) 


DHS showed that wall roughness is inefficient in stimulating small wavelength G dr tier 
vortices and so here we shall be concerned solely with the case when the forced vortex 
wavelength is comparable with the boundary layer thickness. For all the work undertaken 
it transpires that the factor q(k ) is merely a multiplicative constant so henceforth we 
shall take q(k) = 1 which corresponds to a delta-function shaped hump in the spanwise 
direction. Corresponding results for other spanwise distributions may be deduced easily. 
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The remainder of this paper is divided as follows. In §2 we determine solutions 
appropriate to the case when F(x) varies on a relatively fast lengthscale. DHS 
demonstrated that, at least in theory, the asymptotic forms of these solutions can be use 
as the (necessarily unique) initial conditions for marching the full linearised equations 
(1 6 „ b) in the streamwise direction *. These initial conditions are studied in §3 and the 
subsequent numerical solutions addressed in §4. We conclude with a brief discussion of the 

implications of this study. 

2. Vortices induced by short scale roughness elements 

We now examine solutions of (1.6a, 6) when the forcing function F(x) varies on a fast, 
say 0(e), lengthscale. Suppose that the forcing begins at X — x and define 

X = (x — x)/e and F(x) = f(X). (2-D 

DHS pointed out that we may take * = | without loss of generality so that the original 
lengthscale L has been fixed in terms of the distance from the leading edge to the position 
where the forcing begins. In the neighbourhood of * = $ then convective and vertical 
diffusion effects are comparable in a region of depth 0( £ i) and the wall forcing implies 
that u will be 0(1) there. If we define 

(. = y/e’ 


then for £ = 0(1) the disturbance quantities develop as 

(u,u,n.,p)= (n„(X,f),e-U„(JT,0.£-'“o(^0. £ ' S w( X ’0) + "’ (2 ' 2<l) 

whilst / , , \ toon 

(u,v)= (e 5 A«,e 3 pf s ) +•••• * 2 ' 26 ' 

If the governing equations (1.6a, 6) are Laplace-transformed with respect to X and 
solutions of the forms (2.2a) sought then DHS demonstrated that the transformed functions 

U 0 ,v o are given by 




(As) 5 f 


+ AAi 


t) 0 (.s,O = -3A Mf(s) 


(As) 5 f] /•__*_{/ uo(s,q)Ai ((As) 3 <z) dq 
0 | Ai ( ( As ) 5 p)] 

i ' 

(sA)-5 |Ai' [(As)^] - Ai'(O) j + £ J f Ai [(A.s) 3 q 


(2.3a) 


(2.3 b) 
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where a denotes the transform variable and Ai is the usual Airy function. Furthermore, 
by considering these expressions as £ -► oo, it was shown that the flow in this wall layer 
induces a motion in the zone y = 0(1) and there 


U — Uo£ 3 +..., v = Vo£ »+..., 
where U 0 ,V 0 , the transforms of U 0 and V 0 , are given by 

tT 3 F(s)u' . i 
U o = —~(sX) 3 Ai'(0 )?*(!/,*) + ..., 

with q* satisfying the stationary Rayleigh problem 


(2.4a) 


V 0 = -3F(s)sl \ 3 Ai'(0 )q*(y,k) + 

(2.46) 


“ ( d l ~ k 2 ) <?* - «V = o, 9 *(o) 

DHS restricted their attention to the case when F* 
leads to a solution in which 

u ° ~ jffio (</**) , wo ~ 3^r«o (fAV*) ■ (2.4.J) 

This in turn suggests a form of similarity solution as appropriate initial conditions for 
a numerical solution of the system (1.6a, b) to 0(1) distances downstream of the delta- 
function shaped hump. However, DHS were unable to obtain similarity solutions of the 
type assumed in (2 Ad) and so they proposed an alternative form which had the unusual 
feature of predicting that the disturbed flow remains confined to a thin wall layer even at 
0(1) distances from the roughness. We demonstrate that this conclusion is incorrect (see 

§3) but first we examine the inversion of expressions (2.3a, b) for an illustrative selection 
of roughness functions F*(X). 

Since the transformed quantities u 0 , v 0 involve Airy functions of argument proportional 
to ss it is clear that in order to obtained the inverted forms u 0 , u 0 it is essential to cut 
the complex s-plane in an appropriate manner. It is a standard result of transform theory 
that if Q(s) is the Laplace transform of any suitably well behaved function Q(X) then 

i yc+ioo 

Q{X) = 27 , 

where the line of integration passes to the right of any singularities in Q(s). If we cut 
the 5-plane along the arg(s) = tt ray then we can invert Q(s) by integrating around the 
classical ‘key-hole’ shaped contour sketched in Figure 1 and which is composed of a large 
semi-circular arc of radius R and linear portions £i,£ 2 ,£ 3 ,£ 4 and C 5 . The first two of 
the straight line contours are located just above and below the branch cut and a small 
circle C, radius 6, which has centre 5 = 0. For all physically realistic roughness shapes it is 


= 1, ?*(oo) = 0. (2.4c) 

= S(X) and showed that this choice 
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straightforward to verify that, assuming the semi-circular arc is chosen so as to avoid any 
singularities in Q(s), then as R - oo the integrals over the straight lines C 3 ,C* and the 

semi-circle — > 0 and the inverse of Q(s) follows once 6 — > 0. 

We now demonstrate how contributions from C 2 , C and any poles can play 
contrasting roles as the forcing function F'(X) is varied. Interest is restricted to the 

case X » oc for this limit indicates the appropriate initial conditions for a numerical 

solution of equations (1.6a, b ) in order to obtain the vortex characteristics for I = 0(1). 

2.1. The delta function F*(X) = S(X). 

Now F(s) = 1 and so u 0 and €- 0 have no poles in the -s-plane. In addition the 
contribution from C — 0 as its radius shrinks to zero so that the dominant part of 
the solution must arise from the integrals along £1 and C 2 . If we define 

-sX = 8\ ( = X^/\K ( 2 - 5 ) 


9 A '5 at 

u ° = 0 - v T N ' 
2mX 3 


(2.6a) 


N( t,) = B 4 exp(-0 3 )[e^ Ai' (eric*) - c Ai' (fye » ) - Ai'(O) 

+8 2 V l°° (e - ^ Ai (dte*) - e 1 ^ Ai (dte~*)) dt\d6. 
Jr) 


Elementary analysis demonstrates that 

N "" + lifN'" + InN" = o, 


( 2 . 66 ) 


JV(0) = JV'(O) = 0, 1V"(0) = ^fh(J), JV"( 0 ) = ^r(l), (2.6c) 

and where ' denotes a derivative with respect with to r,. It is possible to express N in 
terms of the Whittaker function W K%V defined in Abramowitz & Stegun (1965) so that 

N" = -^Ai(0)r(— D^expl-^r/ 3 ) W 2 ^ (|i/ 3 ) • ( 2 - 7 ) 

9 > 2 

A study of the inversion of (2.3a) leads to the solution 

-- - sM-itfS ( 7Ai<0) + ® 

+ |i)V" - jjv' + ^I) 2 iv - pexp (-Ip 5 ) / (| + ^t* - ) AT(t) ex P (jf )dt , 


(2.8) 
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aj.d the solutions (2.6a), (2.8) together give the large X solution (2.2a). Although they 
are not required in this study, the relevant forms of the spanwise velocity and pressure 
isturbance components, w 0 and p 0 , may now be retrieved from (1.6c, d). 

2.2. The Heaviside function F*(X ) = H(X). 

In this case F( S ) = 1/s and the inversion of u„(s,f) no longer has a negligible 
contribution from the small circle C as 6 -a 0. If the substitutions (2.5) are made it 
IS a straightforward but lengthy task to deduce that as X -> oo 


vo = -■ 


93A-3 

4tt 2 xI 


r (-|)r(|)iv(r/), 


where 


N"" + l^N'" + |r/7V" = 0 : TV( 0) = TV'( 0) = 0, TV — ► 1 

or, equivalently, 

3V3\~h n 


(2.9a) 


as t) oo, (2.96) 


v 0 = -- 


2nXi 

In turn, inversion of (2.3a) leads to 
9 


l (j 


(2.9c) 


u 0 


r (| ){Joo exp (~» p ^ 1 v °W ex p(h 3 ) d <i rfp-i} r&cp(-lt 3 )dt 

r Mhfr .. . i 00 (2-w) 

/ ex p(-^ 3 ) / 


+ / exp(-iy J )| / c 0 (9)exp(I 9 3 )d 9 jdp. 

2.3. A distributed roughness element. 

In order to illustrate the effect of a distributed as opposed to an isolated roughness 
patch we now consider a forcing function given by 


F*(X) = sin X. 


(2.10) 


In this case F(s) = 1/(1 + ,s 2 ) and this introduces simple poles into the Laplace inversion 
integrals required in order to deduce u 0 and v 0 . For this particular choice of forcing 
unction F(s) -* 1 as « -4 0 and thus the contributions to u 0 (X,{), v 0 {X,{) arising from 

ie parts C u C 2 and C of the contour sketched in figure 1 are identical to those relevant 
to the delta function shaped hump problem considered above. 

Inclusion of the terms arising from poles in F(s) leads to the complete solution 


V Cl = (1F± 


3A: 




(2.11a] 
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as 


X — ► oo where N(rj) is as defined in (2.6) and M(X,() satisfies 


d 2 M < x +*f. 


d? 


= — e 


3 Ai 


with 


DM 2 i 

M(x, o) = o, -^r(*’°) = ^J cosX - 


(2.116) 


(2.11c) 


Furthermore, if we denote the solution (2.8) of the delta function problem by u„(X, 0 


then 


u 0 (X,O = {u 0 {X,£) + ...) + |A 3 (Q + (X,0 + Q-(X,£)) 

+ sm ( Ai(Aieif ° eiX _ Ai(X ’ e '* ()e ~ ,x ) • 


(2.12a) 


with 


Q ± (X,t) = e*« X + 



Ai(A»e ± '» () 
3Ai(0)As 




l -r i-2- 
“ 7 ? ^ ' 3 


t)dt + A 3 e 


Ai'(0) Ai(A 



(£Ai(A ±e±%t)dt) 
|Ai(A3e ±! ^ u)| 


-du j. 
(2.126) 


Although the periodic forcing (2.10) is very specific it does illustrate much of the nature 
of the Laplace transform solutions that arise from a wide variety of distributed functions. 
If the forcing has a zero mean value then F(s) tends to air 0(1) constant as « - 0 and so 
the £j , £ 2 and C parts of the contour for the Laplace inversion gives a solution proportional 
to (2.6a) for e 0 and (2.8) for « 0 . On the other hand if F(s) has a simple pole at the origin 
(which corresponds to a distributed forcing function with non-zero mean) then we recover 
a multiple of the Heaviside type solution presented in §2.2. Each pole m F(s), other than 
at the origin 3 = 0, leads to terms like M(X,() in (2.11) and Q ± (X,() in (2.12a). These 
functions are periodic in X and for a general periodic forcing function these terms in the 
solutions for v 0 and u 0 need to be replaced by infinite sums akin to Fourier senes. 

The solutions (2.11a), (2.12a) indicate that as X -+ 00 the influence of the contributions 
arising from any poles of F(s) is confined to a thin region adjacent to the wall whereas 
those from the delta-function part of the solution expand into a region where ( = 0(X 3 ). 
Thus the solution parts divide into forms governed by two distinct structures and this 
point is explored further in the coming section. 
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3. The evolution of 0(1) wavelength vortices further downstream 

We now show how the solutions of the preceding section may be used to deduce unique 
initial conditions of a form suitable for the marching of equations (2.6a, b ) to 0(1) values 
of x. As we allow X — * oo the similarity variable implied by the results of §2 is 


II 

X ~ x — x^ 

(3.1) 

and for ( — 0(1) the appropriate expansions for the disturbance velocity 
and v are 

components u 

As 

u = -~-MQ + • • • , 

X 

v = TT v o(C) + ---, 

X 3 

(3.2 a, b) 

in the case of F*(X) — 6(X). The functions u 0 and vq satisfy 


d 2u ° i 1 duo , d 4 v o 

d ( 2 + 3 ^ d( + ^ U ° “ V °' d ( 4 

+ i X 3v ° + i „ 
3 ^ dC 3 dC ~ ’ 

(3.2 c, d) 

which need to be solved subject to 



dv o 

u o = v 0 = -j— =0, on £ = 0; 

dv o 

«o, ► 0 as ( — *• oc 

(3.3) 


DHS observed that equation (3.2d) may be solved for d 2 v 0 /dC, 2 in terms of Whittaker 
functions. This equation is linear and the constant of proportionality in the solution 
comes directly from the Laplace solution (2.6a) so that 



Now DHS argued that the homogeneous form of equation (3.2c) has the eigensolution 
U ° H = Cexp (-|C 3 ) an d that any other linearly independent solution of (3.2c) behaves 
like u 0 ~ C 3 as C oo. In turn, since u o //(0) = 0 it was stated that the inhomogeneous 
form of (3.2c) can only be solved if a finite multiple of the algebraically decaying solution 
is retained. In that case the solution cannot be matched with the flow in the core region 
where y = 0(1). It was therefore concluded that this similarity solution cannot be the 
required form for x small. Instead, DHS proceeded to identify another similarity solution 
in which to leading order u 0 ~ x^UotC) and u 0 ~ i~iv 0 (C) and these forms of solutions 
led to the conclusion that the disturbance is confined to the region where £ = 0(1). 

However the above argument is flawed as a particular solution of (3.2c) may be found 
which excludes the algebraically decaying homogeneous solution whilst at the same time 

vanishing at ( = 0. (That this is possible was kindly pointed out to us by Professor S.N. 
Brown.) Then 


u 0 =DCexp(— ±( 3 )- 


+ K 


- fsC^o 


+ Cexp(-!C 3 ) f (| + ^< 3 )a 0 (i)exp(I< 3 )d<, 

Jo 

(3.5) 
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is a suitable solution of (3.2a) where D is some constant whose value is as yet unknown. 
Its value could be determined by examining further terms in the similarity expansions 
(3.2a, b) but it is much more convenient to deduce the required value directly from the 
Laplace transform solution (2.7). This gives 

A 4 •! / /%\ \ 

(3.6) 


D = 


1 / , , 3 Ai (0) 

— 7=— 7 Ai(0) H — T7ni 
2y/3n V Al (°) 


) r(|). 


Solutions (3.4-6) describe how the vortex develops at small 0(1) distances downstream of a 
delta-shaped hump and within 0(x 5 ) distances of the wall. They therefore provide suitable 
starting conditions for a numerical solution of equations (2.6a, b) to greater distances 
downstream. 

For the case F*{X) = H(X) the Laplace solutions (2.9a, c) suggest that for 0 < x <C 1 


with 


Now 


d?u o 

w 


u = U 0 (C) + • • 


— X 3 Vo (C) + • • • 7 

(3.7) 

1 >2 du 0 

+ ^lc =V0 ' 

d 4 v 0 

d C 4 

+ i c *' pva + <-( <Pvo -o. 

(3.8) 

3\/3A"5 [< 



(3.9) 

1 

to 

V* ) 



in order to match with solution (2.9c). There is a subtle difference in the solution of the u 0 
equation (3.8a) in contrast to that of equation (3.2a). In the latter case it was seen that of 
the two homogeneous solutions of the equation one decays algebraically whilst the other 
decays exponentially and also vanishes at C = 0. Hence the unknown D was contained m 
solution (3.5) and needed to be deduced by using extra information. However (3.8a) has 
the homogeneous solutions u„ = 1 and u 0 = exp (-£*»)<« so this equation has a unique 

solution when the boundary conditions «„(0) = 1, «o - 0 as C - °° are imposed (with 
the first condition coming from (1.6e)). Thence 


u 0 


= A y C exp(-|t 3 )df + f ex v(-^{J 0 u o(g)exp(|g 3 )dqjdf, (3.10) 

where the constant A is chosen so that uo(0) 1. 

In the two cases described above we have derived solutions of equations (1.6a, b ) valid 
for 0 < x < 1 and y = 0(x *). To deduce the requisite forms for y = 0(1) we let the 

similarity variable ( —* oo and then 

A t«' 1 


u = ^=-«*(v,*) + . 

x^u 


fl+/? 


X 


q*(y, k ) + 


(3.11) 


I f° r 


where the function q*(y,k) has been defined by (2.4) and the constant 8 
F*(X) = 6(X) and (3 = for the Heaviside forcing. In this way we have solutions 
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valid for small distances downstream of the roughness element and these may be used to 
initiate the numerical solutions of (2.6a, 6). 

Finally we discuss the case of the distributed roughness element. As before we shall 
concentrate on the particular case F(X) = sin X where X = i/s although the analysis is 
easily extended to other functions. The solution for X « 1 is composed of two parts. The 
firs of theses just the delta-function solution given in §3.1 above whereas the second 
contains the fast streamwise variation arising from the form of the forcing function. Then 
where the similarity variable C = 0(1) the solutions for u and a develop according to 


(3.12a) 


( Ah£ \ / / ~\ \ 5 \ 

T “•«) + -) + (( T ) + 

" ~ (iT*^ 0 + ■) + (s _ *Vb(jf,c) + ...), 


(3.126) 


where u 0 ,u 0 are as given in (3.2-4) and 


+ Vo=0, ~( - o 

dx de[dx ~ °- 


(3.13) 


These equations cannot be solved subject to the necessary conditions that U 0 = 

, Sm and ~ 0 on y = TIlen 111 an inner layer where Y = y/ £ h =0(1) the 
fast part of the solution (3.12) develops as 


+ ut(X,Y ), 


+*-*v 0 t (*,n 


(3.14) 


Vo = P+e iX + P_e~ ix , 


P± - ^ [_, 


xy? A + = 


A _ 

0 &y 2 ’ 


(3.15a) 


J roo 

Y Ai (-' i ‘ ± ^OoK-A-je^(Ai' ( Aie±?r ) -Ai'(0))] (3.156) 


(3.15c) 


subject to U 0 (X, 0) = sin X and u\ - 0 as Y - oo. It is straightforward to verify that 
as i — > oo so 

V„U3A-Ui'(0)co.,(Jf-Ix) and u) 3A~* Ai'(O) cos(X + !*)/<;, (3.16) 

which shows that the functions U 0 (X,() and V„(X,Q defined in (3.12) are given by 


T y _ 3A“5 Ai'(O) 


s(X + Itt), Vo = 3A 5 Ai'(0)cos(A: - I*). 


(3.17) 
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Finally we have that in the y = 0(1) region the ‘slow' and ‘fast’ parts of the solution 
develop according to 


u ^£^ 4 *( y ,i) + + dtf„(x,y) + . 

X 5 

v~ -W(v.*0 + + <r S V' 0 (X,y) + ... 


(3.18) 


where 

C 0 = 3*-* Ai'(0),-(y,*)cos(X - *.)■ = 3A-JS' Ai'(0),-(y.<c)cos(X + W/« 

and o*(y k) has been defined by the Rayleigh equation (2.4). 

aIK We consequently have the result that for the distributed forcing problem the solut.cn 

structure in the vicinity of the start of the forcing is composed of three d.st.nc reg.on, 

The component of the solution associated with the fast streamw.se 

Rayleigh equation for y - 0(1). is given by the of equa 

c = 0(1) and is governed by viscous-type balances where y - 0{e )• 

contrasts to that relevant to the isolated forcing case in which this last reg.on adjacent 

y = ForO(l)Wluesof X we can seek a solution of (1.6) corresponding to a distributed wall 
roughness with non-zero mean value. If y = 0(1) then 

u = us(*,v) + eW(*,A,y) + • • • > (3.19) 

v = v s (x,y)+e~iv F (x,X,y) + 

„ , „ rp periodic in X with zero mean and are, at leading 

where the fast functions Uf and f 1 , - f /i a\ anc \ 

order in x, proportional to «•(,. *). If the expansions (3.19) are subst.tut d *(W - d 

we average will, respect to the fast variable X we see that „ s and u s also sat.sfy (1.6 . 

We note here the fact that u F and up are driven directly by the wall orcing anc J 

to the forcing in wall layers of thickness 0(d). Moreover, we stress that the fact 

tL Gbrtler munber is 0(1) meal, that no growth or decay of the vortex .s poss.b e on 

the short 0(e) scale so we monitor the evolution of the vortex just ... terms of j* * 

Since the small x forms of these functions are given by the leading order arms n ( • - 

. . j. j?*( y - '. civ. X forcing is, at leading order m 

it follows that the stability problem for the F (X) - sin A g , 

• 1 t; r „1 to that for the delta function case. In a similar way, for distnbuted peno 

functions which provoke the small slow x response appropriate to a Heaviside toct.o.Ohe 
“ problem defined in terms of the slow response of the vortex , .dent.cal to t„ 

, . «,/ x) _ H ( X) It follows that, so long as we define instability in en ■ 

for the case F (A ) - H{X ). p c&n me results from the cases 

of the appropriate downstream scale, . . b r „ , vort i cef; 

F * (X ) = 6(X) and H(X) in order to deduce the stability proper les 

induced by a general wall forcing shape. 
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4. The numerical work 

The similarity solutions derived in the preceding section yield the unique vortex forms 
valid at small distances downstream of the forcing function. In order to investigate the 
development of vortices at 0(1) downstream distances it is necessary to solve equations 
(1.6a, b) numerically. Since these equations are parabolic in nature a standard marching 
scheme was implemented in order to follow the flow from x=l + Ax where Ax is a small 
specified distance from the start of the forcing. The details of the numerical technique 
closely follow those given by Hall (1983) and the reader is referred to that article for a full 
description. Certain modifications were found to be required for the current problem- in 
particular equations (1.6a, b ) were finite-differenced in y but in order to resolve the initial 
vortex form which is concentrated close to y = 0 it was necessary to replace the uniform 

grid employed by Hall (1983) by a variable spaced mesh with nodes tightly packed near 
y ~ 0 but relatively sparsely spaced as y oo. 

The solutions presented in this paper were all verified by making the usual checks 
that the results were independent of grid spacing and choice of Ax. As a measure of the 
disturbance energy we monitored the value of 


-jf 


(a 2 + v 2 + w 2 )dy, 


and defined the local growth rate 0 = E~'dE/dx. The position of neutral stability can 
be taken to be the location at which 0 = 0 and this position does depend to some extent 

on the choice of the particular flow property used to define 0. (This point is discussed in 
detail by Hall (1990).) 

For a selection of values of k with fixed values of G we calculated the position at 
which the induced vortex structure begins to grow. The local wavenumber k x and local 
G or tier number G x were evaluated at that point and used to generate a neutral curve in 
the (k x , G x ) plane. It is worthwhile to emphasise that since the forcing was applied at the 
same point (x = ±) for each calculation, the Gdrtler number G is a measure of how close 
the forcing is to the unstable region. 

We first report on the calculations for the delta- function forcing considered in §2.1. 
In this case the leading order similarity solution valid a small distance downstream of the 
forcing at x = I is given by (3.2a, b), (3.4)-(3.6), (3.11) and these were used to initiate 
the computations at x = \ + Ax. However we found that these leading order solutions 
are insufficiently accurate for our purposes for when these were used our computed neutral 
curves were not independent of grid spacing or the choice of Ax. Therefore it was deemed 
necessary to develop a further term in the similarity solution (3.2). This is a routine task 
and it was found that these two leading terms in the small (x - I) solutions were sufficient 
to ensure adequate convergence. 
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Our results for the delta-function problem are summarised in Figure 2 where we 
show neutral curves for a variety of values of G. We observe that the lowest critical 
Gortler number occurs when G is about 8 and then takes a value of about 15.5. This is 
somewhat greater than that found by DHS (who quoted a lowest critical value of 12) whose 
computations were started with the incorrect similarity forms. This therefore suggests 
that localised wall roughness in its most dangerous form produces growth m the vortex 
structure when the local Gdrtler number is roughly 15.5. In Figure 3 we give further 
details of the envelope of the neutral curves of Figure 2. We see in Figure 3 that initially 
G c decreases with G up to the point where G « 8 after which G c increases with G. The 
critical wavenumber k c is a monotonic decreasing function of G and at the lowest critical 

G c , k c « 0.92. 

In Figure 2 we have also indicated the asymptotic form of the neutral stability curve 
which is valid for large k c . This asymptotic expression may be deduced from the workings 
of Hall (1982) and we see that the first two terms of this asymptotic expansion gives 
surprisingly good agreement with our computed values even at quite modest values of k c . 

We turn now to study the problem of Heaviside forcing as developed in §2.2. In this 
case we again initiated our calculations using the first two terms in the requisite similarity 
solution forms in order to achieve satisfactory convergence. The dependences of critical 
Gortler number G c and local vortex wavenumber k c are shown as functions of G in Figure 
4. We note the radically different behaviours of these forms when compared with the 
results for the isolated delta-function forcing problem. In particular we note that both k c 
and G c are monotonic functions of G and, in particular, G c —■ * 0 as G — > 0. This means 
that the critical Gortler number can be reduced to an arbitrarily small value by choosing 
the origin where the forcing begins increasingly close to the leading edge. Thus it appears 
that when the forcing commences near the leading edge only a very small amount of wall 

curvature is required in order to cause vortex growth. 

As was pointed out in the preceding section we may deduce the vortex response due to 
a general distributed wall forcing by examining the stability characteristics of the Gortler 
modes for the delta- function and Heaviside wall problems. When the forcing function 
has zero mean value then the appropriate stability characteristics close to the start of the 
forcing are reminiscent of the results for the delta-function case whereas for a non-zero 
mean valued wall forcing it is the results summarised in figure 4 which are pertinent. 


5. Conclusion 

In this work we have described how isolated and distributed roughness elements can 
stimulate the growth of linearised Gortler vortices. In particular, Laplace transform 
techniques are required in order to deduce how vortex modes evolve close to the forcing. 
Integration of the full governing equations allows the development of vortices to be 
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monitored at 0(1) distances downstream of the commencement of the forcing. This 
numerical integration requires suitable initial vortex profiles to be specified and these 
(necessarily unique) forms may be found in the form of similarity solutions which may be 
deiived directly from the Laplace solutions. Two specific forcing functions have been 
considered in our work. First, a delta-function was taken as a model of an isolated 
roughness element and the stability characteristics of induced vortices computed. This 
calculation was previously performed by DHS but faulty initial conditions meant that 
theii conclusions regarding the receptivity of 0(1) wavelength vortices were unreliable. 
Our improved workings have led us to conclude that isolated elements can trigger unstable 
vortices at local Gbrtler number greater than about 15.5 but for G x less than this value 
the flow is not susceptible to vortex instabilities. 

When the isolated element is replaced by a distributed roughness patch we have shown 
that the stability characteristics of induced modes are radically altered. Of importance, 
we have seen that as the Gortler number decreases to zero so the critical local wavenumber 
k c grows and the critical local Gbrtler number G c -> 0. This suggests that distributed wall 
forcing can be extremely important in generating Gbrtler modes since irrespective of the 
particular value of G x unstable modes can exist. As was pointed out by DHS, for free- 
stream disturbances the critical G c is roughly 6 so that free-stream disturbances are clearly 
less dangerous in provoking Gbrtler vortices than are distributed wall elements. However 
these free— stream disturbances are more significant than isolated roughnesses and, for a 
general roughness profile composed of both distributed and short-scale variations we have 
shown that it is the former type which tends to dictate the receptivity properties of the 
flow. DHS suggested that the localised wall roughness mechanism would be relevant only 
in an experimental facility with remarkably tiny free— stream disturbances. Our work here 
allows us to extend this conclusion and speculate that a distributed wall element provides 
the easiest route to the generation of vortex modes. 

Unfortunately there are few experiments carried out on the Gbrtler vortex problem; in 
fact we are unaware of any experiments which address the receptivity problem for Gbrtler 
vortices. The fundamental experimental difficulty with the problem is that results appear 
not to be reproducible from day to day with the same apparatus. It is likely that this is 
caused by the sensitive dependence on initial conditions which the vortices exhibit, see Hall 
(1983) for an explanation of why this is the case. However if wall roughness is sufficiently 
large then it is probable that it will lead directly to the linear growth of vortices, the 
analysis we have given here can then be used to predict the response of the flow. Note here 
that it is a simple matter to extend the analysis to describe three-dimensional objects, see 
DHS, so that controlled experiments could be carried out in order to see the response of 
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the flow to different wall disturbances. Thus we argue that the present calculation gives 
the experimentalist, for the first time, a chance to obtain agreement between theory and 
experiment in the G or tier problem. 
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Figure 1 . The key-hole contour used to evaluate the Laplace inversions of § 2 . The complex 
s- plane is cut along the negative real axis and the contour consists of two large arcs of radii 
R, linear portions lines £ \ . £2 , £3, £4 and £5 as shown. The line £1 and £2 are located 
above and below the branch cut and a small circle C of radius 8. Integrals along the line 
parallel to the imaginary axis passing through the point s = c + 0/' may be deduced by 
allowing R — > 00 , 8 —> 0 and then applying Cauchy’s residue theorem. 
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G x 



Figure 2. The neutral curves corresponding to a delta-function shaped wall forcing and 
an initial disturbance constructed from the first two terms of the similarity solution which 
starts as in (3.2). The curves correspond to the values G = 5., 5.7, 6.5, 7.5, 10., 12. 5 t 15., 20. 
The da-shed curves are the one term and two term forms of the appropriate large 
wavenumber right-hand neutral branch (derived from Hall 1982). 







Figure 4 a) The critical local Giktler number G, and b) the critical local vortex 
wavenumber k c as functions of the Gortler number G appropriate to a Heaviskle-function 


wall forcing. 
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